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An overview of the calculation of the two-loop helicity amplitudes for scattering of two gluons into two photons 
is presented. These matrix elements enter into the recent improved calculation of the QCD background to Higgs 
boson decay into a pair of photons, which is the preferred search mode at the LHC for the case of a light Higgs 
boson. 


Recent years have seen an enormous improve¬ 
ment in our ability to calculate two-loop ampli¬ 
tudes with more than a single kinematic invari¬ 
ant. The initial calculations of this type were 
for four-point scattering in maximally supersym¬ 
metric theories [IB and special helicity configura¬ 
tions in QCD [ In this talk we will focus on 
the two-loop amplitudes for gluon fusion into two 
photons [ ||] . This calculation is among the more 
general two-loop processes 

that have become doable, thanks to the develop¬ 
ment of rather general two-loop integration tech¬ 
niques [|i- Further summaries of these develop¬ 
ments may be found in a number of talks at this 
conference [Ill- 

Gluon fusion into two photons is phenomeno¬ 
logically interesting because the preferred mode 
for discovering a light Higgs boson {Mh < 140 
GeV) at the LHC is through its decay into two 
photons. There are large backgrounds to this de¬ 
cay coming from radiation from either partons 
or hadrons [0 [T^ . Although the gluon fusion 
contributions to the background are formally of 
higher order in the perturbative expansion they 
are greatly enhanced by the large gluon distribu¬ 
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tion at small x at the LHC. 

The amplitudes described here were calculated 
in the helicity formalism, using the ’t Hooft- 
Veltman dimensional regularization scheme [ 0 - 
Helicity methods have a long history, having been 
used rather productively at tree and one-loop lev¬ 
els (see e.g. refs. [ 00 )- At two loops the hrst 
of the 2 —> 2 amplitudes were also evaluated using 
helicity states [|J- In the gg —> 77 process un¬ 
der consideration here, the tree amplitudes vanish 
and the one-loop amplitudes give the leading or¬ 
der contributions. Thus the next-to-leading order 
contributions to gluon fusion require an interfer¬ 
ence of two-loop amplitudes with one-loop am¬ 
plitudes. Instead of evaluating this interference 
directly, we again make use of helicity. 

To generate the loop momentum integrals, we 
did not use Feynman diagrams, but instead used 
a unitarity-based technique [ 0,1^. This tech¬ 
nique exploits a duality between loop and phase 
space integrals and has already been employed 
in a number of two-loop calculations [ p], p|, 0. 
This duality was also used very recently to cal¬ 
culate the exact NNLO contribution to the total 
cross section for Higgs boson production at the 

LHC [|g. 

In performing the loop momentum integrals 
some minor extensions of the general reduction 
algorithms developed for four-point massless in¬ 
tegrals [ 1^ were needed. Any four-dimensional 
polarization vector can be expanded in terms of 
the three independent momenta in the problem 
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plus a dual vector. This means that dot products 
of polarization vectors with loop momenta can be 
re-expressed in terms of dot products of loop mo¬ 
menta with external momenta and also the dual 
vector. Integrals where no dual vector appears 
can be reduced using the previously constructed 
reduction algorithms. This leaves integrals con¬ 
taining dual vectors dotted into loop momenta to 
be evaluated; as discussed in ref. [ ^ such inte¬ 
grals tend to be simpler to evaluate. An alterna¬ 
tive method for dealing with helicity at two loops 
was recently given in ref. [ . 

We take the quarks in the loops to be massless 
since the mass of the Higgs boson, and therefore 
the energy scale of the experiments for which this 
calculation is relevant, is well above the mass of 
all quarks other than the top. Moreover, the top 
quark can be ignored since the aforementioned 
energy scale is well below the 2mt = 350 GeV 
threshold. 

The renormalized gg —> 77 amplitude may be 
expanded as 
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where = 3 in QCD and Nf is the number of 
light flavors. In this formula the infrared diver¬ 
gences are encoded as poles in the dimensional 
regularization parameter e = (4 — D)/2. 

The one-loop amplitudes were given in ref. [ ^ 
through their relation to the one-loop four-gluon 
amplitudes. We can write 




where Qi is the electric charge of the zth quark 
and is equal to the sum of non-cyclic per¬ 

mutations of the fermionic contributions to the 
four-gluon primitive amplitudes [ ^. 

It is convenient to extract overall spinor phases 
from each helicity amplitude, 
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where A4g^l^-y^ is the Lth loop contribution, 
as{fM^) is the MS running QCD coupling and a is 
the QED fine structure constant. 

We use Catani’s ‘Magic Formula’ [ |2^ for two- 
loop infrared divergences to organize our results. 
This formula was obtained for general QCD am¬ 
plitudes, but with minor modifications it also ap¬ 
plies to mixed QED and QCD amplitudes. In 
cases where the tree-level amplitude vanishes (as 
happens for gg 77 ), Catani’s formula reduces 
to, 

Ai(2) =/(i) (e)A4*^^^ (1) 
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where (e) contains the infrared singularities, 

and Afd)®'’ is a finite remainder. In our case the 
color factors can only be proportional to ^ gQ 
/d)(e) 

is relatively simple: 
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trough 0{e^) the amplitudes simplify greatly 
the m[^\ , , reduce to, 
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where we are using an “all-outgoing” convention 
for the momentum (p^) and helicity (A^) labeling. 
The Mandelstam variables are s = (pi -I- ^ 2)^1 
t = {Pi+ P 4 )^, and u= {pi+ pa)^. 

We consider both QCD corrections with inter¬ 
nal gluon lines and QED corrections with internal 
photons. For the QCD corrections, the depen¬ 
dence of the finite remainder in eq. (|^ on quark 
charges, N, Nf and the renormalization scale p, 
may be extracted as. 
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The two-loop renormalized QED corrections are 
a little simpler, since in this case the amplitudes 
are free of infrared divergences. 


x5a, 


pSL 

A2A3A4'^ Ai A 2 A 3 A 4 


(3) 


We quote our results in the physical s-channel 
(s > 0; t, u < 0). In order to reduce the size of 
the expressions we define 
t u 

x = -, y=-, X = ln(-a;), E = In(-p), 

5 S 

X = X -\- ZTT, Y = Y -\- ZTT, 

S = X2-k7r2, T = E2 + 7r2, 

Z±=X±Y, Z = (X - E)2-h 7^^ 

TT^ 

= Li„(-x) ± Cn, B = Li 2 (-a;) - —, 
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Ct{x,y) = Li„(-a;) ±Li„(-p). 


The explicit forms for the appearing 

in eq. (^) are 
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+ ^(2x‘^ - 54x - 277/2) (-1 + -f ) . 
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Similarly, the subleading color contributions in 
eqs. (H) and i) are, 
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FfL_+is,t,u) = F5+_+(s,77,t). 

The reliability of these results was ensured by per¬ 
forming a series of checks described in ref. [ p|. 

A companion talk in these proceedings [ 
describes the application of the amplitudes pre¬ 
sented here to obtain [ Q an improved prediction 
for the QCD background to Higgs production at 
the LHC, when the Higgs decays into two pho¬ 
tons. 
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